SCIENTIFIC 

REPpRTS 



n5 




OPEN 



SUBJECT AREAS: 

CLIMATE SCIENCES 

CLIMATE CHANGE 

CRYOSPHERIC SCIENCE 

CLIMATE AND EARTH SYSTEM 
MODELLING 

Received 
26 February 201 3 

Accepted 
13 June 2013 

Published 
2 July 2013 



Correspondence and 
requests for materials 
should be addressed to 
S.H. (stephan. 
harrison@exeter.ac. 

uk) 



Rapid thinning of the late Pleistocene 
Patagonian Ice Sheet followed migration 
of the Southern Westerlies 

J. Boex 1 , C. Fogwill 1 ' 2 , S. Harrison 1 , N. F. Glasser 3 , A. Hein 5 , C. Schnabel 4 & S. Xu 6 

College of Life and Environmental Sciences, University of Exeter, Exeter, EX44RJ, UK, 2 Climate Change Research Centre, University 
of New South Wales, Sydney, New South Wales, Australia, institute of Geography and Earth Science, Aberystwyth University, 
Aberystwyth, Wales, SY23 3DB, 4 NERC Cosmogenic Isotope Analysis Facility, SUERC, East Kilbride, G75 OQF, UK, 5 School of 
GeoSciences, University of Edinburgh, Edinburgh, ScotlandUK, 6 Scottish Universities Environmental Research Centre (SUERC), 
Scottish Enterprise Technology. Park, East Kilbride, UK G75 OQF. 

Here we present the first reconstruction of vertical ice-sheet profile changes from any of the Southern 
Hemisphere's mid-latitude Pleistocene ice sheets. We use cosmogenic radio-nuclide (CRN) exposure 
analysis to record the decay of the former Patagonian Ice Sheet (PIS) from the Last Glacial Maximum (LGM) 
and into the late glacial. Our samples, from mountains along an east-west transect to the east of the present 
North Patagonian Icefield (NPI), serve as 'dipsticks' that allow us to reconstruct past changes in ice-sheet 
thickness, and demonstrates that the former PIS remained extensive and close to its LGM extent in this 
region until —19.0 ka. After this time rapid ice-sheet thinning, initiated at —18.1 ka, saw ice at or near its 
present dimension by 15.5 ka. We argue this rapid thinning was triggered by a combination of the rapid 
southward migration of the precipitation bearing Southern Hemisphere (SH) westerlies and regional 
warming. 

Reconstructing the timing and rate of decay of the Earth's ice sheets since the Last Glacial Maximum (LGM) 
allows us to assess the mechanisms of climate change globally during a major climate transition 1,2 . This 
provides crucial information for testing climate models, and helps us better predict the response of modern 
ice sheets to future climate change. Of particular value are reconstructions of the rate and timing of glacier 
changes regionally and globally during the last glacial to interglacial transition, as these provide insights into 
global climatic hemispheric teleconnections, the mechanisms by which external climate forcing are transmitted 
through the Earth's climate system 3 . 

Establishing if the onset of deglaciation was synchronous globally allows detailed assessment of the relative 
roles of atmospheric and oceanic circulations as drivers of climate change. An atmospheric driver is suggested 
from high resolution synchrony between the Earth's climatic systems of the northern and southern hemispheres 
and implied from the comparison of glacial and palaeoecological records at many mid latitude sites which closely 
follow changes recorded in the Greenland ice core records 4,5 . Alternatively, an oceanic driver of climate signalling 
at this time is suggested from comparison with Antarctic ice core records, and mid to high latitude sites in 
Southern Patagonia, suggesting that during the glacial-interglacial transition there may have been asynchronous 
behaviour in the climate systems of the two hemispheres 6 " 8 . 

Resolution of the debate is important because of the light it sheds on the mechanisms of climate change. If 
glacier and palaecological changes on millennial and sub -millennial scales are in phase in both hemispheres, then 
this points to the dominance of atmospheric circulation mechanisms of climate change. If, on the other hand, the 
fluctuations are out of phase, then it suggests the operation of an oceanic bipolar seesaw linked to the thermoha- 
line circulation that causes warm conditions in the north to coincide with cold oceanic conditions in the south and 
vice versa 8 . Other, regional factors may also play a role, including the suggestion that regional warming was 
caused by C0 2 release from the Southern Ocean triggered by insolation-driven Antarctic sea ice retreat 9 . Testing 
these hypotheses will improve our understanding of the sensitivity of these ice sheets to past and future climate 
forcing, and is critical to understanding the global drivers of climate change 3 . In order to do this, however, a more 
complete record of the behaviour of mid-latitude Southern Hemisphere ice sheets through time will be needed 
than has been achieved so far. In working towards this goal we present here the first three-dimensional recon- 
struction of a major outlet of one of the former Southern Hemisphere (SH) mid-latitude ice sheets. Our focus is 
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Figure 1 | Location map; Map A, South America and the location of the study area also showing the existing North and South Patagonian Icefields 
(dark grey) together with the mapped LGM ice sheet (light grey) 11 . Map B is adapted from 17 and shows the modelled ice thickness based on mapped outer 
limits in the Lago Buenos Aires Lobe. The white box shows the location of the study area and transect that forms the focus of this study (see Figure 2). 



the Pleistocene PIS during the LGM and understanding its evolution 
provides insight into the role of the SH climate system during degla- 
ciation. 

The Patagonian Icefields consist of two separate ice masses, the 
NPI and South Patagonian Icefield (SPI). They are the largest tem- 
perate ice masses on Earth and their outlet glaciers are some of the 
most dynamic. The icefields are nourished by precipitation from the 
Southern Hemisphere westerlies and their dynamics are controlled 
by seasonal variations in the westerlies and associated ocean cur- 
rents 410 . Geomorphological and chronological evidence shows that 
during the late Quaternary the Patagonian Icefields coalesced to form 
the PIS over the Southern Andes 11 (Figure 1). Over the past eighty 
years researchers have mapped and dated the limits of the ice sheet 
providing firm constraints on its lateral extent during multiple glacial 
phases throughout the Quaternary 11 " 14 . This makes it the ideal loca- 
tion to assess the forcing mechanisms of climate driven glacial 
change in the Southern Hemisphere. However whilst past work 
has established the extent and timing of ice sheet advances these 
lateral constraints do not allow us to reconstruct changes in the 
vertical extent and volume of the ice sheet 14 " 16 . This latter informa- 
tion is crucial to help us understand the nature of climate forcing and 
ice sheet response, and the contribution of the ice sheet to past sea 
level rise. 

Geomorphological evidence indicates that fast-flowing outlet gla- 
ciers occupied both the Lago Buenos Aires valley and Lago 
Pueyrredon basin; two major ice discharge routes of the PIS. In the 
Lago Pueyrredon basin two moraine systems, the Rio Blanco and 
Hatcher moraines, have been dated by CRN methods to the LGM 
and MIS 8 (260 ka), respectively 1315 . The Rio Blanco moraines were 
deposited between 29 ka and 21 ka and mark the horizontal LGM 
extent in this area 131516 . The age for the Hatcher moraines was deter- 
mined from outwash cobbles to 260 ka 15 . Fast-flowing glaciers 
played a key role in effectively drawing down the main central ice- 
mass with outlet lobes characterised by low surface gradients and 
inferred low basal shear traction associated with basal sliding. High 
spatial resolution three-dimensional ice sheet models of central 
Patagonia 17 at the LGM show a 'highly dynamic, low angled ice sheet' 
which was drained by large ice streams to the east and west and with a 
mean ice thickness of —1130 m. The modelled ice sheet was con- 
sidered 'thin' when compared to previous regional -scale modelling 
reconstructions of the PIS at the LGM 18 . The modelling result 17 
agrees with geomorphic interpretations showing fast-flowing outlet 



glaciers drained the eastern margins of the former ice sheet 11 
(Figure 1). 

Whilst the horizontal ice limits are reasonably well established for 
large parts of the LGM PIS, there are no reconstructed ice limits to 
constraint its vertical extent or ice surface profile. This has hindered 
reconstructions of past ice sheet volume and rate of decay and our 
ability to test the impact of climate on ice volume through time, 
critical to understanding ice sheet dynamic processes. To derive 
detailed chronological constraints on late Quaternary ice surface 
elevation changes we targeted mountains thought to have protruded 
through the ice sheet above the Pueyrredon Lobe 17 . Our study site 
forms a west-east transect extending from Cerro Tamango (1722 m; 
47°10' 06.87S; 72°34'30"W) in the west, through Cerro Oportus 
(2076 m; 47°07'14"S; 72°07'53"W) to Sierra Colorado (1537 m; 
47°22'05"S; 71°37'23"W) (Figure 2). The mountains separate the 
Lago Pueyrredon basin to the south from the Chacabuco Valley to 
the north. We used geomorphological analysis and CRN exposure 
analysis on erratic boulders, bedrock and moraine boulders deposited 
on the mountain flanks to reconstruct the west-east ice-surface pro- 
files. The northern flanks of Cerro Tamango and Cerro Oportus were 
sampled and these locations record upper ice limits in the Chacabuco 
Valley. In light of recent CRN calibration rate studies in Patagonia we 
apply the production rate for the isotopes 10 Be and 26 A1 derived from 
New Zealand that overlap at 1 sigma with an independently derived 
production rate from Lago Argentino, Patagonia 16,19 . We present the 
data across the profile in Figure 2, together with our interpretation of 
the former ice sheet surface. The full data set is available in Table SI 
in the Supplementary Information. 

Results 

The uppermost constraint on ice sheet thickness comes from Sierra 
Colorado where two distinct lateral moraines mark the former ice 
surface (Figure 2). 10 Be and 26 Al CRN exposure ages of 113 to 177 ka 
were obtained from 3 erratic boulders on the uppermost moraine at 
1368 m altitude. The exposure ages are similar to those obtained 
from boulders on the Hatcher moraines further east 13 . Based upon 
the analysis of paired 10 Be and 26 A1 we suggest that the true age of this 
moraine is in excess of 177 ka. The result supports previous mapping 
and indicates the moraine represents an expansion of the PIS that 
took place before the last interglacial 15 . The valley floor to the south 
sits at around 150 m elevation and thus the PIS was around 1,200 m 
thick here at this time. 
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Figure 2 | Digital elevation model (SRTM) of the Lago Pueyrredon basin showing the study area, transect and reconstructed ice surfaces. The 

white lines in the east (upper panel) show previously mapped outer limits, the Hatcher and Rio Blanco moraines 15 . Pink dots show the location of 
cosmogenic exposure age samples, from mountains Cerro Tamango in the west to Cerro Oportus and Sierra Colorado in the east. The white dotted line 
shows the location of the transect shown in the cross section. The lower panel shows the transect generated from SRTM data with xl 5 vertical exaggeration 
(see Figure 1 for location). Sampled boulders are shown with exposure ages (circles show 10 Be, triangles show both 10 Be and 26 A1 exposure ages). 



The age of the lower moraine identified on Sierra Colorado, 100 m 
below the uppermost ice- limit gives exposure ages ranging from 
24.5-28.9 ka. Based on geomorphological evidence we suggest the 
oldest sample reflects the most likely age of this moraine. This agrees 
with the CRN estimates of 28.7 ka for the outermost of the Rio 
Blanco moraines, 190 km from the centre of the contemporary 
NPI and 20 km to the east of Sierra Colorado 1516 (Figure 2). These 
data demonstrate that both the Rio Blanco moraine and the lower 
moraine on Sierra Colorado date to the LGM as previously sug- 
gested 1314 . Our results indicate that the PIS reached a maximum 
elevation of 1100 m at the LGM in the vicinity of Sierra Colorado, 
120 km from the centre of the NPI. This estimate agrees well with 
high-resolution ice-sheet modelling reconstructions 17 and is sup- 
ported by erratic CRN exposure ages of 17.9 ± 0.7 ka and 18.9 ± 
2.6 ka from the summit of Cerro Oportus, and this demonstrates the 
summit remained covered by ice at the LGM. The profile of the ice 
sheet at this time can be tied to a lateral constraint, the Columna 
Moraine System 13 (Figure 2). Based on an exposure age from an 
erratic on the Columna Moraine system and the exposure ages from 



the summits of Cerro Oportus we infer that the ice sheet extended 
horizontally over 120 km from the centre of the modern NPI at 
around 20.0 ± 2.4 ka. 

Our dating indicates that the summit of Cerro Oportus at 1895 m 
altitude was exposed at around —19.0 ka, with a stepped lowering of 
the ice surface to 1300 m altitude at —18.1 ka. After this time we 
record rapid exposure of the full altitudinal profile of Cerro Tamango 
and Oportus, suggesting rapid thinning after 18 ka, with over 
1000 m of vertical thinning of the former ice sheet within approxi- 
mately 1000 years. This interpretation of rapid deglaciation is well 
constrained by three combined 10 Be/ 26 Al analyses, which are intern- 
ally consistent (Figure 2). After this period of rapid ice sheet draw- 
down, there was a brief and limited readvance, recorded at the Maria 
Elena Moraine at around 17 ka, and this demonstrates that ice was 
present in the Chacabuco Valley at this time. 

Discussion 

Our constraints on the LGM and lateglacial deglaciation record 
shows the PIS response to climate forcing through this important 
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Figure 3 | Deglacial records illustrating the rapid thinning of the 
Patagonian Ice Sheet at 47.1 °S against warming and rising atmospheric 
C0 2 in Antarctica, upwelling in the Southern Ocean and NGRIP ice core 
record 22 , (a) Time-series of ice sheet volume and area during deglaciation 
from the optimum LGM extent at 23,500 through to 1 1,000 years driven by 
ELA re-scaled from the Vostok temperature reconstruction 22 Periods of 
warming in Antarctica are highlighted in red, while the Antarctic Cold 
Reversal (ACR) is highlighted in blue and the Younger Dryas (YD) is in 
grey, (b) Alkenone-based SST reconstruction from core MD07-3128 
(53 °S) 28 . (c) Accepted exposure ages for ice sheet stages, marking the 
changing volume and extent of the ice sheet in the Lago Pueyrredon basin 
and the Chacabuco Valley calculated with the independently derived NZ 



production rate 1619 , (d) Opal flux from ocean cores TN057-13PC (51°S, 
4°E) and NBP9802-6 (62°S, 169°W) as a proxy for upwelling in the 
Southern Ocean 22 . Opal flux and February SST are plotted on the published 
age model for TN057-13PC 32 . (e) February SST estimated by applying the 
modern analog technique 30 ' 31 to diatom species assemblages extracted from 
TN057-13PC 33 ' 34 . (f) Ice isotope chronologies for the EPICA Dome C 
(EDC) and EPICA Draunning Maud Land (EDML) core (light and dark 
blue respectively), and North Greenland Icecore Project (NGRIP) core 
(orange) 35 . 

climate transition, with initial retreat after 29.0 ka and the ice sheet 
remaining close to its LGM dimension until — 19 ka (see Figure 3 A- 
E). Data from Cerro Oportus and Cerro Tamango indicate stepped 
thinning between 19 ka and 18.1 ka, followed by rapid thinning 
throughout the altitudinal profile until a brief still stand or advance 
at around 16.9 ka. Evidence from a series of former ice-dammed 
lakes within this limit (Figure 3D) suggests that ice had withdrawn 
rapidly to within 10-15 km of its present extent by 15.6 ka 13 There is 
no evidence here of a substantial readvance in the Lago Pueyrredon 
basin related to either the Antarctic Cold Reversal (ACR) or The 
Younger Dryas (YD). Subsequent expansion of the ice sheet during 
the late glacial served only to dam westwards drainage to the Pacific, 
leading to the establishment of a series of lakes filling the basins to the 
east of the present day NPI 11,13 . 

This direct record of large-scale rapid thinning of the ice sheet 
between 19.0 and 15.6 ka clearly highlights the sensitivity of the PIS 
at this latitude to changing climate (Figure 3a, c). Our observations 
have two important implications. Firstly, they support high resolu- 
tion ice-sheet models 17 with respect to the rate of ice sheet decay 
through the Last Glacial- Interglacial Transition (Figure 3F), and 
emphasise the key role fast flowing outlet glaciers played in effec- 
tively drawing down the main central ice mass 17 . Secondly, even with 
the uncertainties in the data the timing of this thinning is coincident 
with marked changes in both the SH oceanic and atmospheric sys- 
tems. We can suggest several hypotheses that can account for this 
pattern of ice sheet thinning. In particular thinning coincides with 
marked warming of SH mid to high latitudes, that has been linked to 
changes in ocean thermohaline circulation, specifically warming of 
the southern ocean —the bipolar seesaw— thought to be triggered 
either by spring insolation changes combined with variations in 
greenhouse gas fluxes 20 or by ocean circulation changes due to 
increased melt water from the Northern Hemisphere ice sheets at 
ca. 19 ka 2 ' 20 ' 23 ' 27 . Although SST changes between MIS 4-2 played a 
significant role in driving the evolution of the PIS, its behaviour 
lagged SST. The reasons for this are unclear but probably relate to 
changes in ice dynamics 21 . However, concurrent to this regional 
warming there is also physical evidence 22 of a rapid southward lat- 
itudinal shift of the precipitation bearing Southern Hemisphere 
Westerlies (SHW). This is particularly evident in cores in the SE 
Pacific at this time 23,24 (Figure 3A-D). 

Regional evidence for this marked atmospheric shift comes from 
marine records south of the Antarctic Polar Front (53.2°-61.9°S) off 
the southern Chilean margin that record increased opal accumula- 
tion and Alkenone-based reconstructions linked to enhanced upwel- 
ling during the late glacial (Figure 3b and d). It is suggested that this 
enhanced upwelling was triggered as the core of the SH westerlies 
migrated rapidly southwards from their LGM position. Whilst 
debate remains over the exact timing marine records suggest that 
southern migration was initiated at around —17 ka 24 , and the SH 
westerlies shifted rapidly south and stabilised at —15.5 ka, by which 
time they may well have reached 62 °S in the SE Pacific 1,25 , 
(Figure 3d) 23 . A southward shift in the SH westerlies would have 
warmed the Southern Ocean and Antarctica 22 by allowing degassing 
of C0 2 to the atmosphere and providing a positive feedback to initial 
warming 26,27 . This poleward shift would have also significantly 
reduced precipitation in the vicinity of the current NPI, and we 
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propose that this reduction in precipitation in conjunction with 
regional warming amplified the rapid ice sheet decay that we record 
(Figure 3e,f). This suggestion is further supported by climate mod- 
elling studies that predict significant increases in precipitation at the 
LGM in NW Patagonia 24 . However, recent work 28 has found that 
oceanic changes associated with the bi-polar seesaw was sufficient 
to explain glacier recession and ELA rise in southern Patagonia dur- 
ing late glacial times, and they argue against the shift in the Southern 
Westerlies as a trigger for deglaciation. This may not be the case in 
central Patagonia however. Glasser et al. 29 dated Younger Dryas age 
moraines at the mouths of several valleys around the NPI, suggesting 
this argues against a bipolar seesaw operating at these latitudes dur- 
ing this time. Clearly, this issue is not yet resolved 30 . 

Our findings highlight the sensitivity of central Patagonia to cli- 
mate change which during Pleistocene times may have been com- 
municated via regional warming and shifts in the latitude of the SH 
westerlies. We argue that the latitudinal variation in the core of the 
precipitation bearing SH westerlies is probably key to defining the 
extent and magnitude of glaciation at these latitudes, although the 
pattern to the south may have different causes. Our results have 
important implications for glaciers and hydrological systems in 
Patagonia if, as predicted, westerly airflow continues to migrate pole- 
wards as a result of a warming climate 27 . Our approach using moun- 
tains to record the three -dimensions of the former ice mass makes 
important steps towards aligning the geomorphology of Patagonia 
with ice sheet and climate modelling based studies, and provides 
critical new insights on the sensitivity of this region to climate 
change 31 . 

Methods 

In order to reconstruct the three-dimensional evolution of the PIS in the study area we 
used geomorphological analysis to ground-truth the landform mapping from satellite 
imagery, maps from the Institute Geographical Militar of Chile, Landsat imagery and 
aerial photographs. The chronology of ice sheet thinning was obtained using glacial 
erratics and glacially eroded bedrock which were sampled from key landforms and 
altitudinal profiles for CRN exposure analysis using 10 Be and 26 Al (see SOM and 
Figure 2). 

Samples for CRN were used to constrain the age of moraines or ice sheet trimlines 
and also to reconstruct the timing of ice sheet thinning and mountain top exposure. 
Samples were taken from stable boulders and bedrock surfaces and around half a 
kilogram of whole rock was required for analysis. Samples were processed at the 
laboratories at the NERC CIAF at SUERC and at the University of Exeter. Crushing 
was followed by mineral separation and chemical etching to produce clean quartz. 
The quartz was then dissolved in order to chemically extract and separate 10 Be, and for 
seven samples 26 A1. The extracted isotopes were then measured by AMS at SUERC. 
Results for 10 Be and 26 Al concentrations were converted to exposure ages using the 
CRONUS-Earth online calculator, version 2.2 using the NZ Macaulay landslide, NZ 
calibration data set (Supplementary Table S2). The Dunai time varying model is used 
for the data presented here to allow direct inter- comparison with previous studies (see 
SOM for details of this). No correction was made for erosion, snow cover or isostatic 
uplift in this study and therefore the exposure ages presented are minimum ages. 
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